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The LHC observed so far three deviations from the Standard Model (SM) predictions in flavour 
observables: LHCb reported anomalies in —>■ K*pT and R{K) — B ^ Kp^p~/B —>■ Ke^e~ 
while CMS found an excess in h —>■ pr. We show, for the first time, how these deviations from the 
SM can be explained within a single well-motivated model: a two-Higgs-doublet model with gauged 
Lfi — Lt symmetry. We find that, despite the constraints from r — > ppp and Ba-Bs mixing, one can 
explain h —>■ pr, B —>■ K*p'^p~ and R{K) simultaneously, obtaining interesting correlations among 
the observables. 


I. INTRODUCTION 

So far, the LHC completed the SM by discovering 
the last missing piece, the Brout-Englert-Higgs parti¬ 
cle [BE]. Furthermore, no significant direct evidence 
for physics beyond the SM has been found, i.e. no new 
particles were discovered. However, the LHC did ob¬ 
serve three ’hints’ for new physics (NP) in the flavor sec¬ 
tor, which are sensitive to virtual effects of new parti¬ 
cles and can be used as guidelines towards specific NP 
models: h —>• pr, B —>• K*p^p~, and R{K) = H ^ 
Kp^p~/B Ke^e~. It is therefore interesting to ex¬ 
amine if a specific NP model can explain these three 
anomalies simultaneously, predicting correlations among 
them. 

LHCb reported deviations from the SM predictions [3l 
E] (mainly in an angular observable called Pg 0) in 
B —>• K*p'^p~ |5] with a significance of 2-3 cr depending 
on the assumptions of hadronic uncertainties [THg. This 
discrepancy can be explained in a model independent ap¬ 
proach by rather large contributions to the Wilson coeffi¬ 
cient Cg [Mia, i e. an operator {s^aPLb){pj°'p), which 
can be achieved in models with an additional heavy neu¬ 
tral Z' gauge boson [MU]. Furthermore, LHCb [T6| re¬ 
cently found indications for the violation of lepton flavour 
universality in 

which disagrees from the theoretically rather clean SM 
prediction P™ = 1.0003 ± 0.0001 [IB by 2.6 a. A possi¬ 
ble explanation is again a NP contributing to Cg ^ involv¬ 
ing muons, but not electrons [T8U20] . Interestingly, the 
value for Cg required to explain R{K) is of the same 
order as the one required by P — K*p'^p~ 0 EB- 
In Ref. [TB], a model with gauged muon minus tauon 
number — Lr) was proposed in order to explain the 
B —>■ K* p'^ p~ anomaly. 

Concerning Higgs decays, CMS recently measured a 
lepton-flavour violating (LFV) channel [32] 

Br[h ^ pr] = (0.89t°;4°) % , (2) 

which disagrees from the SM (where this decay is forbid¬ 
den) by about 2.4 cr. Such LFV SM Higgs couplings are 


induced by a single operator up to dim-6 and Br[/i —>■ pr] 
can easily be up to 10% taking into account this op¬ 
erator only [23ll28j . However, it is in general difficult 
to get dominant contributions to this operator in a UV 
complete model, as for example in models with vector¬ 
like leptons [33] ■ Therefore, among the several attempts 
to explain this h ^ pr observation [3nH34| . most of 
them are relying on models with extended Higgs sec¬ 
tors. One solution employs a two-Higgs-doublet model 
(2HDM) with gauged — Lr [35] - 
The abelian symmetry 17(1)is interesting in gen¬ 
eral: not only is this an anomaly-free global symmetry 
within the SM [53H55] . it is also a good zeroth-order ap¬ 
proximation for neutrino mixing with a quasi-degenerate 
mass spectrum, predicting a maximal atmospheric and 
vanishing reactor neutrino mixing angle |391141] . Break¬ 
ing — Lr is mandatory for a realistic neutrino sector, 
and such a breaking can also induce charged LFV pro¬ 
cesses, such as T —>• 3/r [33ES] and h ^ pr [3H] - 

Supplementing the model of Ref. |35j with the in¬ 
duced Z' quark couplings of Ref. m can resolve all three 
anomalies from above. Interestingly, the semileptonic B 
decays imply lower limit on g' jMz', which allows us to 
set a lower limit on t — >■ ppp, depending on /i — >■ pr. 


II. THE MODEL 

Our model under consideration is a 2HDM with a 
gauged symmetry [3S]- The L^ — Lr symme¬ 

try with the gauge coupling g' is broken spontaneously 
by the vacuum expectation value (VFV) of a scalar $ 
with = 1, leading to the Z’ mass 

mz' = V2g'{<^) = g'vis >, (3) 

and Majorana masses for the right-handed neutrinos^. 

Two Higgs doublets are introduced which break the 
electroweak symmetry: dti with (3^%^ = —2 and '^2 


^ Active neutrino masses are generated via seesaw with close-to- 
maximal atmospheric mixing and quasi-degenerate masses ESj. 
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with = 0. Therefore, ^2 gives masses to quarks 

and leptons while d'l couples only off-diagonally to r/r: 

— Q— QfY^j^ 2 di + h.c.. 

Here Q (£) is the left-handed quark (lepton) doublet, u 
(e) is the right-handed up-quark (charged-lepton) and 
d the right-handed down quark while i and / label the 
three generations. The scalar potential is the one of a 
{7(l)-invariant 2HDM [33] with additional couplings to 
the SM-singlet $, which most importantly generates the 
doublet-mixing term 


that induces a small vacuum expectation value for 
d>i [33]. We define tan/3 = (5'2)/('I'i) and a is the usual 
mixing angle between the neutral CP-even components 
of 'I'l and '1^2 (see for example |31j)- We neglect the 
additional mixing of the CP-even scalars with Re[$]. 

Quarks and gauge bosons have standard type-I 2HDM 
couplings to the scalars. The only deviations are in the 
lepton sector: while the Yukawa couplings of 'I '2 

are forced to be diagonal due to the — Lt symmetry, 
Cr/j gives rise to an off-diagonal entry in the lepton mass 
matrix: 

/j/e sin /3 0 \ 

® yMsin/3 0 • (5) 

\ 0 ^TfiCOsP yrsin/dj 

It is this T-/r entry that leads to the LFV couplings of h 
and Z' of interest to this letter. The lepton mass basis 
is obtained by simple rotations of (/r_R,T/j) and (/rL,ri) 
with the angles Or and 9l, respectively: 


sin Or ~ 


y2r 


-Ctm cos/3, 


tan 9l _m^ 
tan Or nir 


( 6 ) 


The angle di, is automatically small and will be neglected 
in the following.^ A non-vanishing angle Or not only 
gives rise to the LFV decay h ^ fir due to the coupling 

~ sin(0fl) cos{0R)fPRnh = F(1 fP r/t/i , 

V cos(p) sin(p) ^ 

(7) 

in the Lagrangian, but also leads to off-diagonal Z' cou¬ 
plings to right-handed leptons 


g'K (m. 't) 


/cos 20R 
l^sin 29 R 


sin20R \ 
- cos 20R) 


YPr 



( 8 ) 


while the left-handed couplings are to a good approxima¬ 
tion flavour conserving. In order to explain the observed 
anomalies in the B meson decays, a coupling of the Z' to 
quarks is required as well, not inherently part of — L^ 


^ Choosing ^2 would essentially exchange 0^ 

[35], with little impact on our study. 


models (aside from the kinetic Z-Z' mixing, which is as- 
sumed to be small). Following Ref. [TS], we introduce 
heavy vector-like quarks, i.e. Qr = {Ur, Dr), D^, 
and their chiral partners Qr = {Ur, Dr), Dr, Ur, with 
vector-like mass terms 

mqQRQR + "rnDDhUR + muU rUr + h.c., (9) 

and Lfj, — charges -1-1 (i.e. QrD_r^ = = —1), 

coupling them to the Z' boson. Yukawa-like couplings 
involving the heavy vector-quarks, the light chiral quarks 
and 

3 _ _ 

d> ^ (bRY^PRd, + UrY^PrUj) (10) 

3 _ _ 

+ ^ (^DRY^PRd, + UrY^Pru,) + h.c. 

3=1 

then induce couplings of the SM quarks to the Z' once $ 
acquires its VEV. Thus, integrating out the heavy vector¬ 
like quarks gives rise to effective Z'dtdj couplings |3H1I3H] 
of the form 


5 ' {da^PRd.z'^rf^ + d.^^PRd,z'^Tf^), ( 11 ) 

with hermitian matrices F^/' that are related to the vec- 

‘'J 

tor quark masses and Yukawa couplings 

as follows [15] : 


pdi? ^ _ 


2m|, 


{Y,W*), 


2ml 


(yQvQ*), (12) 


which holds in the approximation |r®p^| <C 1.^ 


III. FLAVOUR OBSERVABLES 

We will now recall the necessary formula in the region 
of interest (i.e. small Or) considering only the processes 
giving to most relevant bounds on our model, i.e. Bg-Bg 
mixing, neutrino trident production and r —3/i. 


A. h ^ rt 

The branching ratio for h ^ rt reads 

BT[h^ rt] - ^ ( 13 ) 

where Fsm — 4.1 MeV is the decay width in the SM for a 
125 GeV Higgs [37] and F(l^ is defined in Eq. 0. Com¬ 
paring this to Eq. ([^ one sees that both sin6*R ^ 0 
and cos (a — /3) ^ 0 are required to explain the CMS 
excess [35] . 


^ Compared to Ref. m, the vector-like quarks also have Yukawa 
couplings to the — Lr-charged scalar doublet This 

induces a small additional mass mixing among the heavy quarks, 
and also a coupling to h suppressed by cos(q;—/3). We assume 
these couplings to be small to avoid large contributions to 
gg ^ h and h 77 . 
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tan(/3) = 10 (yellow) and tan(/J) = 50 (blue) 




FIG. 1: Left: Allowed regions in the cos(a — /3)-sin(6_R) plane. The blue (light blue) region corresponds to the la (2a) 
region of the CMS measurement of /i —>■ /rr for tan j3 = 50; yellow stands for tan /? = 10. The (dashed) red contours mark 
deviations of > rr by 10% compared to the SM for tan/3 = 50 (10). The vertical green lines illustrate the naive LHC limit 
I cos(a — P)\ ^ 0.4, horizontal lines denote the 90% C.L. limit on r —>■ 3/t via Z' exchange. 

Right: Allowed regions in the T^s-mz'/g' plane from B —> K*fj,^fi~ and R(K) (yellow) and Bs mixing (blue). For Bs mixing 
(light) blue corresponds to (rriQ = 15mz>/g') mq = mz' jg'■ The horizontal lines denote the lower bounds on mz'lg' from 
r —>• 3yi for sin(0i{) = 0.05, 0.02, 0.005. The gray region is excluded by NTP. 


B. Lepton decays 

While the Higgs contributions to r —y /r/i/r and r —)■ /r 7 
turn out to be very small in most regions of parameter 
space [33 due to the small lepton masses involved, the Z' 
contributions to r —y 3/r can be sizable |42j and restrict 
The branching ratio is given by 


suppressed Cq^). Cg®® is not affected, which natu¬ 
rally generates violations of lepton flavour universality 
in i? —y /B Ke^e~. We find 


'-'9 


/2 1 

9 ^ 1 pdL(fl) 

V2ml,aGFVtbV* 


(15) 


(14) 

which has to be compared to the current upper limit 
of 2.1 X 10“® at 90% C.L. [JS] obtained by Belle. A 
combination with data from BaBar |49j gives an even 
stronger limit of 1.2 x 10“® at 90% C.L. [SD], to be used in 
the following. For small Or, the branching ratio for r —?► 
/i 7 is proportional to the same combination but 

highly suppressed by 2aI f, and hence not as restrictive. 


C. B^ K*g+g,- and B Kfj.+g.-/B Lfe+e" 

Both B —?> and R{K) are sensitive to the 

Wilson coefficients and While in our 

model the contribution to Cio is suppressed by sin(20yj)) 
(or even sin(20L)), the Wilson coefficients and 
with muons are generated (as well as Cg'^ and the Br 


where we set cos(20i^) = 1. As already noted in Ref. [131 
[3T] Cg^ < 0 and C'g^^ = 0 gives a good fit to data. Using 
the global fit of Ref. [8] we see that at (Ict) 2ct level 

- 0.5 (-0.8) > ReC^^ > (-1.6) - 2.0. (16) 

Interestingly, the regions for Cg^ required by R{K) and 
B —y K* lie approximately in the same region. Fur¬ 

thermore, a good fit to the current data does not even 
require Cg^^ [5] , so we neglect it in the following for sim¬ 
plicity. This can be achieved in the limit mo 
resulting in F®*^ » We will also assume our Cg^ 

to be real for simplicity. Note that our model predicts 
the decay B —y Kfir (recently discussed in Refs. [53] ) 
to be suppressed by Bj^ compared to B —?► Kfi^, while 
B —y AT/ie and B —Kre are forbidden. 


D. Bs—Bs mixing 


^ For conventions see Refs. [S[T9]. 


The interactions of Z' and relevant for B —)■ AT/i+zr 
also contribute to Bg-Bg mixing m- For rriR uiq, 
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we get 


Mi2 


ra' 1 + 


~ 1 + 


1 g ™Q 


167r^ 


52 


6477^ 




(17) 


We require the NP contribution to be less than 15% in 
order to satisfy the experimental bounds m- Due to 
the dominance of the vector-quark Q we can express Pfa 
directly in terms of Cg ^ from Eq. (15) and find the upper 
bounds 

<3.2TeV/|C9^^|, mg < 41 TeV/ICg^H ■ (18) 


Combining Eq. (18 1 with Eq. then gives an upper 
bound of mz' Ig' < 4TeV (6.5 TeV) at 1 cr (2cr). 


E. Neutrino trident prodnction 


The most stringent bound on flavour-diagonal Z' cou¬ 
plings to muons arises from neutrino trident production 

(NTP) [T51155] : 


(Tntp 

CTntp 


1 -f (l -f 4sl 


1 -I- (1 -I- 4s^) 


(19) 


Seeing as our region of interest is in the small On regime, 
the NTP bound is basically independent of the angle 
Or. Taking only the CCER data [54], we get roughly 
"mz'!g' ^ 550 GeV at 95% C.L. Compared to r — ggp, 
the trident neutrino bound only dominates for very small 
values of 9r, roughly when Or < 10“^ (see Fig. ^ (right)). 

For mz' > mz, the LHC constraints from the process 
pp —)■ pipiZ' —>■ 4/r (or 3p, plus missing energy) jSS] are 
currently weaker than NTP but will become com¬ 
petitive with higher luminosities [56H58j . 


F. Phenomenological analysis 

Concerning the phenomenological consequences of of 
our model, let us first consider the implications of > 
pLT. In the left plot of Fig. [^we show the regions in the 
cos(a — /3)-sin(0jj) plane which can explain h ^ pir at 
the 1 cr and 2 tr level for different values of tan /3. Mea¬ 
surements of the h couplings to vector bosons require 
I cos(a — /3)| ^ 0.4 |59l |60] while the Higgs effects in 
r —> 3/x and r —>• ^7 are typically negligible [3S]. As 
a side effect, the h ^ pir rate also implies a change in 
the h ^ TT rate, although this is negligible in regions 
with small Or. In addition we show the regions com¬ 
patible with T —>■ 3^ for various values of mz'jg' ■ Note 
that g' ^ 0.3 in order to avoid a Landau pole below the 
Planck scale. In summary, small values of Or can explain 
the CMS h ^ pLT excess for moderate to large values of 
tan/3 for cos(a — j3) ~ 0.1. 

In the right plot of Fig. [^we examine which regions in 
parameter space can account for B —>■ K* pT taking 
into account the constraints from Bg-Bg mixing. Since 


cos(q:-/3) = 0.2 


0.10 


0.08 


0.06 


cc 

05 


c 

'c/3 
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0.00 
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FIG. 2: Allowed regions in the mz'/g'-sm{9ii) plane: the 
horizontal stripes correspond to h ^ pir (la) for tan/3 = 
85, 50, 25 and cos(q;— /3) = 0.2, (light) blue stands for (future) 
r —>■ 3p limits at 90% C.L. The gray regions are excluded by 
NTP or Bg-Bg mixing in combination with the 1 a range for 
C 9 (see Eq. @). 


we focus on the limit Mr —)• 00 (i.e. Cg —)■ 0 ) we find 
that unless Pgl" is rather large, B K*pL^pT can be 
explained without violating bounds from Bg-Bg. Only a 
very small Pgg independent region is excluded by NTP. 
In addition, bounds from t ^ 3pi depending on sin(dii;) 
can be obtained. 

Concerning r —>■ 3//, future sensitivities down to 
Hr [t —)■ 3/t] ~ 10 “® seem feasible El] and will cut deep 
into our parameter space (see Fig. Using the 1 a lim¬ 
its on h ^ piT to fix 9r and Bg mixing with Cg to fix 
mz'!g' - as well as the LHC limit | cos(a — /3)| < 0.4 - 
we can obtain a lower limit on the rate r —>■ 3/t 

Hr [t 3/i] > 3.8 x 10"® (10/ tan pf , (20) 


which implies tan /3 > 18 with current data m and 
tan/3 ^ 61 if branching ratios down to 10“® can be 
probed in the future. This is the main prediction of our 
simultaneous explanation of > /tr, B —)■ K*pC and 
R{K). 

Finally, we remark that a Z-Z' mixing angle Ozz' jUj 
is induced by the VEV of 4/i |35| 


\g'0zz>\ 


giv"^ cos^ /3 
m^Jg'^ 



TeV y 

mz'lg'J 

( 21 ) 


which leads to small shifts in the vector couplings of Z 
to muons and taus 


gvigg, tt) ~ -1/2 -h 2s^ ± g'Ozz’l{glcw) , (22) 

and thus ultimately to lepton non-universality |43j . For 
the values of interest to our study (see Fig. [^, and in 
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the limit mz ^ wz', the shift is automatically small 
enough to satisfy experimental bounds and leads to tiny 
branching ratios Z —>■ /ir below 10“® (for Or < 0.1). 
Note that the couplings to electrons and quarks remain 
unaffected. For mz' S> mz, the p parameter is enhanced 
by gSl 



and is therefore compatible with electroweak precision 
data (p — 1 < 9 X 10“^ at 2cr [51]) for the parameter 
space studied in this letter. 

IV. CONCLUSIONS 

In this letter we showed for the first time that all 
three LHC anomalies in the flavour sector can be ex¬ 
plained within a single well-motivated model: A 2HDM 
with a gauged — Lr symmetry and effective Z'sb cou¬ 
plings induced by heavy vector-like quarks. Except for 
the T-p couplings, the Higgs sector resembles the one 
of a 2HDM of type-I. Therefore, the constraints from 
h decays or LHC searches for A° t+t“ are rather 
weak and h ^ pr can be easily explained in a wide 
parameter space. The model can also account for the 
deviations from the SM in B —>■ K*p~ and naturally 
leads to the right amount of lepton-flavour-universality 
violating effects in R{K). Due to the small values of the 
T-p mixing angle Or, sufficient to account for h ^ pr, 
the Z' contributions to r —>■ 3/i are not in conflict with 


present bounds for large tan /3 in wide rages of parameter 
space. Interestingly, B K*p'^p~ and R{K) combined 
with Bg-Bg put a upper limit on mz' /g' resulting in a 
lower limit on r —>■ 3p if Br[/i —>■ pr] ^ 0: for lower val¬ 
ues of tan /? the current experimental bounds are reached 
and future sensitivities will allow for a more detailed ex¬ 
ploration of the allowed parameter space. The possible 
range for the — Lt breaking scale further implies the 
masses of the Z' and the right-handed neutrinos to be 
at the TeV scale, potentially testable at the LHC with 
interesting additional consequences for LFV observables. 
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